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摘 要 : 等 高 线 曲 率 是 一 个 重要 的 地 形 属 性 ,反应 了 地 形 表 面 在 水 平方 向 的 凹凸 性 ,表达 了 地 表 物 
质 运 动 的 发 散 和 汇合 模式 。 基 于 安 塞 县 县 南 沟 小 流域 的 矢量 等 高 线 数据 和 DEM ,分 别 利 用 圆 拟 合 
算法 ( 相 邻 三 点 法 、 间 隔 三 点 法 和 最 小 二 乘法 ) 和 曲面 拟 合 模型 (下 模型 .S 模型 和 2Z 模型) 提取 等 
高 线 曲率 ,通过 对 实地 地 形 的 对 比分 析 , 结 果 表 明 :(1) 在 矢量 等 高 线 数据 的 计算 结果 中 ,三 点 拟 合 
法 计算 结果 相 比 最 小 二 乘法 结果 差异 大 ,对 等 高 线 曲 率 空 间 格 局 分 布 描述 更 准确 ;(2) 最 小 二 乘法 
计算 的 结果 频数 分 布 集中 程度 最 高 ,两 种 三 点 拟 合 法 计算 结果 频数 曲线 差别 微小 ;(3) 在 栅 格 数字 
高 程 模型 的 计算 结果 中 ,基于 S 模型 计算 结果 在 空间 格局 上 较 卫 模型 和 Z 模型 的 结果 差别 大 , 基 


F OE 模型 的 计算 结果 对 等 高 线 曲 率 描述 效果 更 好 。 结 果 能 准确 说 明 采 用 不 同 算法 计算 等 高 线 曲 


率 的 差别 ,对 在 实际 数字 地 形 分 析 中 有 重要 的 意义 ,可 为 等 高 线 曲率 计算 提供 重要 参考 。 


关键 词 : 


地 面 曲率 是 对 地 形 表 面 上 一 点 扭曲 变化 程度 的 
定量 化 度量 因子 '" 1。 等 高 线 曲 率 是 通过 该 点 的 等 
值 面 (水 平面 ) 与 地 表 交 线 的 曲率 ,等 高 线 曲率 表达 
了 地 表 物 质 运 动 的 汇合 和 发 散 模式 。 与 坡度 、 坡 
向 等 基本 地 形 参数 一 样 , 等 高 线 曲率 在 数学 定义 上 
是 明确 且 唯 一 的 ,具有 确定 的 表达 式 "1。 目 前 ,对 
坡度 . 坡 向 和 流域 网 络 分 析 等 的 研究 比较 深入 ,而 对 
等 高 线 曲 率 等 地 形 曲率 的 研究 较 少 。 基 于 栅 格 数据 
结构 提取 分 析 曲 率 , FLORINSKY "通过 误差 传播 
定律 对 常用 的 四 种 曲率 模型 进行 了 分 析 比 较 , 通 过 
计算 在 中 心 点 处 各 个 方向 的 偏 导数 来 计算 曲率 的 精 
度 , 给 出 了 曲率 误差 和 数字 高 程 模 型 (DEM ) 误差 的 
关系 式 。SCHMIDT 中 通过 实验 操作 在 光滑 的 数学 
曲面 和 实地 DEM 数据 对 3 种 曲率 计算 模型 进行 了 
研究 。 但 FLORINSKY 和 SCHMIDT 在 探讨 等 高 线 
曲率 方面 的 研究 较 少 ,更 缺少 对 矢量 等 高 线 的 研究 。 
而 基于 矢量 数据 结构 提取 分 析 曲 率 , 张 亮 提出 圆 
拟 合 的 方法 ,并 没有 进行 详细 对 比 不 同 拟 合 方式 之 


数字 地 形 分 析 ; DEM; 等 高 线 曲率 ; 圆 拟 合算 法 


设计 新 的 基于 矢量 数据 格式 的 等 高 线 曲 率 计 算 方 
法 ,以 期 为 DEM 分 辩 率 的 确定 提供 科学 依据 ,弥补 
前 人 研究 的 不 足 , 且 对 推进 数字 地 形 分 析 的 研究 具 
有 重要 意义 。 


1 等 高 线 曲 率 提取 算法 


1.1 基于 矢量 数据 的 等 高 线 曲 率 算法 

1.1.1 三 点 所 合法 ”如果 三 点 不 在 一 条 直线 上 则 
可 确定 唯一 的 圆 ,中 间 节 点 处 的 曲率 用 半径 的 倒数 
表示 ;如 果 三 点 共 线 , 则 直接 赋 曲 率 值 0 。 

如 图 1 所 示 为 等 高 线 的 一 部 分 ,假设 7 为 三 点 
(Pii Pi Pi ) 构 成 的 三 角形 外 接 圆 半径 , 则 半径 
的 倒数 1/r 即 可 认为 是 在 点 已 处 的 曲率 HEHE 
海伦 公式 .正弦 定理 .三 角形 面积 公式 即 可 推导 出 等 
高 线 曲率 计算 公式 ; 


(a+b t+c)(a+b-c)(ate-b)(b+c-a) 


abe 


C,=1/r= 


(1) 


间 对 结果 的 影响 。 基 于 上 述 原因 ,本 研究 利用 不 同 
算法 提取 等 高 线 曲率 并 进行 对 比分 析 和 评价 ,讨论 
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式 中 :C, 表示 节点 已 处 的 等 高 线 曲 率 (rad © m~’); 
a, b 和 c 分 别 表示 线段 P,_P,, P,P; ? P,P; 69 
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图 1 三 点 拟 合 法 计算 等 高 线 曲 率 示意 图 


Fig.1 Three-point fitting method for calculating 


contour curvature 


图 2 最 小 二 乘法 计算 等 高 线 曲 率 示意 图 


Fig.2 Least squares method for calculating contour curvature 


长 度 。 

1.1.2 最 小 二 乘法 最 小 二 乘法 是 利用 最 小 化 
误差 的 平方 和 寻找 最 佳 拟 合 函数 匹配 该 组 数据 。 如 
图 2, 样 本 集 (X;,Y;) i= (1,2,3…N) 为 等 高 线 上 一 
组 节点 ,d; 为 样本 集中 点 到 圆心 (4,B) WERK, R A 
拟 合 圆 的 半径 , 取 半 径 倒 数 1/XR 即 为 在 所 选 样本 集 
中 的 等 高 线 曲率 。 根 据 圆 曲 线 方程 .最 小 二 乘法 推 
求 拟 合 圆 半径 ,得 到 利用 最 小 二 乘法 求 得 的 等 高 线 
曲率 公式 为 : 


C, =1/R=2/ Ja +P -4e (2) 
, HD - EG HC -ED 
WR: @ = Cop D- CG’ 


YG +Y) +a% X, + Y, 
z= 一 N s 
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E=NYX+NYV XD RR- ¥ (+ YS) DY, 
G=(N> = SYSY), 
H=NYXYVY,+4NYY - DR t+) YY, 
1.2 基于 格 网 DEM 的 等 高 线 曲率 算法 

曲率 计算 的 实现 主要 包括 对 地 形 曲面 求 一 阶 
导 、 二 阶 导 和 混合 导 , 对 于 DEM 数据 曲率 计算 的 实 
现 , 是 利用 局 部 曲面 拟 合 或 者 对 局 部 窗口 进行 差分 
HARP 

通过 3 x 3 局 部 窗口 如 图 3 中 的 各 个 高 程 点 
(XYZ) G =1,2, n), H T RREAN, 
以 最 小 二 乘法 来 计算 多 项 式 中 的 系数 ,然后 按照 等 
高 线 曲率 计算 公式 


2 一 2p0s + pt 
C, = aa ae (3) 
f af af af af af 
: = ’ = ’ = ? t= > = 
式 中 :p ax q ay r x? ay? s axdy 


基于 格 网 DEM 提取 等 高 线 曲率 ,因为 对 地 形 曲 
面 拟 合 所 采取 的 方法 不 同 ,会 导致 计算 的 5 个 参数 
P, gr, s, t 的 值 有 所 不 同 ,从 而 会 使 最 终 计算 的 等 
高 线 曲率 值 有 所 差别 ,本 文 主要 探讨 3 种 方法 拟 合 
地 形 曲面 而 对 提取 的 等 高 线 曲 率 产 生 的 影响 。3 种 
方法 分 别 为 二 次 曲面 方法 (EVANSI 方法 )、 限 制 
二 次 曲面 方法 (SHARY 2” 方法) 和 不 完全 四 次 曲 
MDE (ZEVENBERGEN! 方法 )。 
1.3 等 高 线 曲 率 的 算法 实现 
1.3.1 基于 矢量 数据 的 等 高 线 曲 率 计 算 本 人 研究 
基于 Net 环 境 下 的 C# 程 序 开发 技术 编写 程序 ,提取 


图 3 3x3 局 部 窗口 示意 图 


Fig.3 3 x3 partial window diagram 
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等 高 线 曲 率 。 具 体 等 高 线 曲 率 提取 流程 如 图 4, 主 
要 步骤 为 :(1) 对 原始 矢量 等 高 线 数据 进行 拓扑 检 
查 和 错误 修复 ,包括 对 等 高 线 交 叉 . 自 相交 REE 
等 方面 的 错误 修复 。( 2) EH ArcMap 中 的 Feature 
Vertices To Points 工具 将 拓扑 修复 后 的 矢量 等 高 线 
图 层 转换 为 点 图 层 , 利 用 Add XY Coordinates 工具 
为 等 高 线 节点 图 层 添加 坐标 。(3 ) 利用 C# 编 写 程 


序 ,以 添加 了 节点 坐标 的 点 图 层 中 的 “ORIG_ID" 字 
段 (存储 了 节点 所 属 等 高 线 的 索引 号 ) 为 基础 读 入 
同一 条 等 高 线 上 节点 的 坐标 数据 ,用 上 文 推导 的 公 
式 (1) 和 (2) 分 别 计算 等 高 线 曲率 并 存储 于 点 图 层 
属性 表 的 curvature 字段 中 。(4) 对 矢量 点 格式 曲 
率 结果 进行 空间 分 析 , 完 成 利用 不 同 圆 拟 合算 法 基 
于 矢量 数据 的 等 高 线 曲率 的 提取 与 分 析 。 


等 高 线 曲 率 .txt 


矢量 等 高 线 数据 | 拓扑 修复 > 矢量 等 高 线 数据 “Ve 等 高 线 节 点 图 层 
| 国保 存 为 txt 格 式 等 高 线 节 点 图 层 Add XY Coordinates 
a 节点 图 层 La 
ArcMap 处 理 分 析 - hit 


图 4 ”基于 矢量 格式 的 等 高 线 曲 率 提取 流程 图 


Fig.4 Flow chart of contour curvature extraction which base on vector data 


1.3.2 基于 格 网 DEM 的 等 高 线 曲率 计算 图 5 
为 基于 格 网 DEM 的 等 高 线 曲率 提取 流程 图 :(1) 在 
ArcMap 中 加 载 原始 DEM 数据 ,利用 Toolbox 中 Ras- 
ter to ASCH 工具 将 原始 DEM 数据 转化 为 文本 格式 
的 DEM 文档 。(2) 利用 C# 编 程 实现 3 种 方法 对 等 
高 线 曲 率 的 计算 ,得 到 文本 格式 的 等 高 线 曲率 文档 。 
(3) 将 文本 格式 的 曲率 在 ArcMap 中 利用 转换 工具 
ASCII to Raster 转换 成 栅 格 格式 曲率 。(4) 对 栅 格 
格式 曲率 结果 进行 空间 分 析 ,完成 利用 不 同 算法 基 
于 DEM 等 高 线 曲 率 的 提取 与 分 析 。 


DEM 数 据 HRastertoAscbjDEM 文 本 文档 |VSs 计 算 曲 率 庆 等 高 线 曲 率 .txt 


ASCII to Raster 
栅 格 格式 等 高 线 = 

Eee A — 
图 5 基于 DEM 的 等 高 线 曲 率 提取 数据 流程 图 
Fig.5 Flow chart of contour curvature extraction 


which base on DEM 
2 试验 区 与 基础 数据 


选择 位 于 陕 北 安 塞 县 县 南 沟 流域 的 两 个 试验 区 
的 矢量 等 高 线 数据 和 利用 ANUDEM 软件 插值 生成 
的 栅 格 DEM 数据 ,如 图 6、7、8.9 所 示 。 该 流域 地 处 
黄土 丘陵 沟 密 区 的 第 二 副 区 、 兴 干 旱 森 林 章 原 地 区 。 
该 地 区 海拔 1 070 ~1 429 m, 地 面相 对 起 伏 260 m, 
地 面 坡度 三 30° 达 一 半 以 上 。 该 流域 沟 密林 立 ,地表 
破碎 复杂 ,水 土 流失 异常 明显 ,适用 于 本 人 研究 对 等 高 


| \ 1246m 1 267 m 
Wd \ 
川 J / 
. i 
1 152 m hil 
í R 
1 257 m 1221 m 
0 100 200 1257 m\\ 
Im 2 \ SSH 


1: 10 000 等 高 线 ， 等 高 距 5 m 
到 6 样 区 1 矢量 等 高 线 数据 


Fig.6 Sample area 1 vector data 


QW AS Hs NAKA 


R EANNA ASN ey Wwe 
0 100 200. VV 827 m=) 


1 : 10 000 等 高 线 ， 等 高 距 5 m 
图 7 样 区 2 矢量 等 高 线 数据 


Fig.7 Sample area 2 vector dData 
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2.5 m 栅 格 尺 十 DEM 


图 8 样 区 1 机 格 DEM 数据 
Fig.8 Sample area 1 DEM 


5 mite R DEM 


图 9 样 区 2 栅 格 DEM 数据 
Fig.9 Sample area9 DEM 


线 曲率 的 提取 分 析 。 
3 结果 与 分 析 


3.1 基于 矢量 等 高 线 曲率 结果 分 析 

实验 数据 1 选 自 图 6 县 南 沟 流 域 的 一 部 分 , 实 
验 数 据 2 选 自 图 7 男 一 块 县 南 沟 小 流域 的 一 部 分 ， 
两 部 分 选取 的 原则 都 是 在 地 形变 化 较为 复杂 的 地 
区 ,从 而 计算 的 等 高 线 曲 率 值 大 小 跨度 范围 相对 
大 ,分 类 效果 较为 明显 ,便于 对 不 同 算法 计算 结 
的 比较 分 析 。 实 验 1,2 计算 的 等 高 线 曲率 结果 如 
图 10,11。 结果 表明 :利用 三 点 拟 合法 计算 的 等 高 
线 曲率 比 利 用 最 小 二 乘法 计算 的 等 高 线 曲率 大 ， 
曲率 在 空间 分 布 上 存在 着 明显 差异 ;同一 种 方法 
(三 点 拟 合法 ) 提 取 等 高 线 曲 率 , 选 点 方式 不 同 ,得 


出 的 等 高 线 曲率 在 空间 分 布 上 差异 不 明显 。 在 等 
高 线 比 较 弯 曲 , 变 化 复杂 的 地 区 曲率 大 ,而 在 等 高 
线 比 较 平 直 , 变 化 较 小 的 地 区 曲率 小 ;在 地 形 特 征 
线 山 谷 线 和 山 浓 线 周围 等 高 线 曲率 较 大 ,而 在 非 
地 形 特征 线 的 地 区 曲率 较 小 ,结果 符合 实地 地 形 
情况 。 

3.2 基于 矢量 等 高 线 曲 率 结 果 统计 分 析 

将 两 个 地 区 的 矢量 点 格式 的 等 高 线 曲 率 值 导 入 
Excel 进行 频率 统计 ,得 到 如 图 12 不 同方 法 拟 合 等 
高 线 曲率 频率 曲线 。 通 过 对 图 12 等 高 线 曲 率 频率 
曲线 分 析 ,结果 表明 :在 3 种 拟 合算 法 计算 等 高 线 曲 
率 中 ,最 小 二 乘法 集中 程度 最 高 , 相 邻 三 点 法 和 间隔 
三 点 法 集中 趋势 差别 微小 ;采用 不 同 拟 合算 法 对 等 
高 线 曲 率 的 计算 ,最 小 二 乘法 计算 结果 与 三 点 法 计 
算 结 果 差 别 较 大 ,结果 值 偏 小 ,曲线 表明 该 流域 曲率 
值 集中 分 布 在 区 间 (0,0.4) 中 。 

对 两 个 样 区 曲率 值 进行 统计 分 析 可 见 ( 表 1)， 
样 区 1 等 高 线 曲率 计算 结果 为 相 邻 三 点 法 与 间隔 三 
点 法 计算 结果 相差 微小 , 样 区 2 等 高 线 曲 率 计 算 结 
果 为 相 邻 三 点 法 与 间隔 三 点 法 计算 结果 咯 有 差别 ; 
在 两 个 样 区 中 ,最 小 二 乘法 计算 的 结果 明显 偏 小 。 
3.3 EF DEM 等 高 线 曲 率 结果 分 析 

基于 不 同 曲面 拟 合 模型 的 样 区 1,2 的 等 高 线 曲 
率 计算 结果 (图 13、 图 14) 可 以 看 出 : 

(1) 基于 S 模 型 计算 出 的 结果 较 基 于 EE 模型 和 
Z 模型 计算 的 结果 在 空间 格局 上 差别 大 ,棱角 分 明 ; 
基于 下 模型 计算 的 结果 较 基 于 Z 模型 计算 的 结 
在 地 形 特征 线 附 近 的 曲率 值 差 别 更 大 ,显示 效果 更 
清晰 。 

(2) 在 等 高 线 曲 率 计算 结果 图 上 能 够 清晰 地 看 
到 有 关 的 地 形 特 征 线 ,并 且 地 形 特 征 线 的 分 布 与 实 
地 情形 相符 合 , 分 布 格局 也 合理 ;等 高 线 曲 率 值 在 地 
形 特征 线 山脊 线 或 山谷 线 周 边 较 大 ,并且 越 靠 近 山 
顶 曲 率 值 越 大 ;曲率 值 的 正 负 反映 了 地 形 的 凹凸 变 
化 ,对 地 形 特 征 线 的 勾 绘 是 正确 的 。 

3.4 基于 DEM 等 高 线 曲率 结果 统计 分 析 

在 ArcMap 中 利用 属性 分 类 工具 对 基于 3 种 曲 
面 拟 合算 法 提取 的 等 高 线 曲率 结果 进行 分 区 统计 分 
析 发 现 :对 于 同一 个 样 区 ,采用 不 同 的 曲面 拟 合算 法 
来 计算 基于 DEM 的 等 高 线 曲率 值 的 区 间 分 布 有 差 
别 ,但 频数 的 分 布 趋势 趋 于 一 致 ,都 在 0 值 邻 域内 集 
中 (图 15 .图 16)。 
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10 实验 1 计算 的 等 高 线 曲 率 


Fig. 10 Contour curvature computed from experiment 1 
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11 实验 2 计算 的 等 高 线 曲 率 


Fig. 11 Contour curvature computed from experiment 2 
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Fig.12 Frequency curve of contour curvature fitted by different methods 
Rl 两 样 区 等 高 线 曲率 统计 表 
Tab.1 Contour curvature statistics for two sample areas 
标准 差 ; pale 
> } M M M 
方法 Min Max Mean (STD) 方法 in ax ean (STD) 
样 区 1 相 邻 三 点 法 0 0.903 28 0.042 836 0.057 136 FE 2 相 邻 三 点 法 0 2.662 044 0.082 330 0.135 489 
间隔 三 点 法 0 0.903 28 0.043663 0.057490 间隔 三 点 法 0 1.461236 0.051528 0.066736 


小 二 乘法 0 0.206 361 0.004040 0.011 317 最 小 二 乘法 0 0.127311 0.000110 0.001 513 
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图 13 样 区 1 等 高 线 曲 率 计算 结果 


Fig. 13 Contour curvature computed from sample area 1 


(a) 二 次 曲面 模型 (了 模型 ) O) 限制 二 次 曲面 模型 ( S 模 型 ) | 


\ 


(Q) 不 完全 四 次 曲面 模型 (Z 模 型 ) 


\ J: N 
mA 、 
一 ‘ 


) 


a 
| E:_797 BA 


加 14 样 区 2 等 高 线 曲率 计算 结 


Fig. 14 Contour curvature computed from sample area 2 
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16 样 区 2 等 高 线 曲率 频数 区 间 分 布 


Fig. 16 Frequency interval distribution of contour curvature from sample area 2 
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进一步 对 样 区 1 和 样 区 2 进行 等 高 线 曲 率 数 值 
的 统计 分 析 可 见 ( 表 2) ,对 于 同一 样 区 ,采用 模型 
计算 的 结果 值 域 和 均值 与 S 模型 结果 值 差别 小 ,但 
与 Z 模 型 计算 的 结果 值 差 别 大 , 且 根 据 标准 差 Z 模 
型 曲率 结果 值 更 加 集中 。 


表 2 两 样 区 等 高 线 曲 率 统计 表 
Tab.2 Contour curvature statistics for two 


sample areas 


方法 Min Max Mean STD 
样 区 1 E -2114 1982 -0.07540 14.760 23 
S 模 型 -2985 1915 -0.08630 21.15734 
Z 模 型 -1703 2166 -0.02685 12.347 80 


样 区 2 下 模型 ” -797 313 
S 模型 ” -791 352 
Z 模 型 ” -298 285 


—17.674 528 103.776 805 
—28.142 857 102.319 521 
-4.503 226 74.582 656 


(a) 样 区 1 基于 矢量 和 DEM 等 高 线 曲 率 结果 符合 图 
rad*(100 m)'00®00 @0@ 1.564 rad-m! A 


高 :1982 0 


0oeoeeee2.662 radm- Ss 
0 40 8 


图 17 同一 地 区 不 同 格式 的 元 数据 等 高 线 曲 率 提取 对 比 


Fig.17 Comparison of contour curvature from different 


formats in the same region 


3.5 基于 矢量 数据 和 基于 DEM 提取 的 等 高 线 曲率 
的 比较 

基于 相 邻 三 点 法 入 模型 计算 的 等 高 线 曲 率 
结果 ,对 矢量 数据 结果 和 栅 格 DEM 结果 进行 比 
较 , 选 取样 区 1 和 样 区 2 地 面相 对 起 伏 较 大 ,地 形 
特征 线 附 近 的 区 域 进行 撩 机 结果 对 比 ,结果 如 图 
17 所 示 。 分 析 图 17 以 及 结合 图 13 .14 可 以 看 出 ， 
基于 矢量 数据 和 基于 DEM 提取 的 等 高 线 曲率 值 
在 空间 格局 上 差别 不 大 , 均 能 够 说 明 等 高 线 的 弯 
曲 变 化 ,在 地 形 特征 线 周边 ,等 高 线 曲率 的 变化 也 
一 致 ;基于 规则 格 网 DEM 计算 出 来 的 结果 更 加 平 
滑 细腻 。 


4 讨论 


本 文 仅 对 基于 不 同 算法 提取 等 高 线 曲率 进行 了 
对 比分 析 , 在 基于 矢量 等 高 线 数据 不 能 区 分 曲率 正 
负 的 情形 下 仍 需 进一步 研究 。 本 文 暂 提 供 解决 思 
路 :以 图 18 所 示 等 高 线 为 例 ,A, ,B, AC, 分 别 表示 
等 高 线 上 三 点 ,A, 和 Ci 两 点 确定 唯一 直线 , 取 线段 
A.C, 中 点 D ,判断 点 B, 处 高 程 值 和 点 D, 处 高 程 
值 大 小 , 即 可 解决 基于 矢量 等 高 线 不 能 区 分 曲率 正 
负 的 情形 ,判断 依据 参照 公式 (4) 。 


Ep <E, WRR +” (Qn) 18 点 B, 和 点 D, ) 


Es > 有 ,曲率 取 ”- ”( 如 图 18 点 B, 和 点 D,) (4) 
E, = 五 ) ,曲率 取 “0” 


式 中 :Es 为 B 点 高 程 值 ;E， 为 D 点 高 程 值 。 


R 


18 示例 等 高 线 


Fig.18 Example contour 


在 后 续 研究 中 将 落脚 点 着 重 放 在 算法 实现 方 
向 ,以 期 能 达到 更 为 理想 的 结果 ,更 好 的 改进 原 有 算 
法 的 不 足 ,为 确定 合适 DEM 的 分 辩 率 提供 重要 
人 参考。 


5 结论 

本 研究 基于 矢量 等 高 线 数据 和 栅 格 DEM ,分 别 
利用 圆 拟 合算 法 和 曲面 拟 合 模型 提取 等 高 线 曲率 并 
进行 对 比分 析 ,通过 设计 新 的 基于 矢量 数据 格式 的 
等 高 线 曲率 计算 方法 ,为 DEM 分 辨 率 的 确定 提供 科 
学 依据 ,弥补 因 基于 等 高 线 派生 的 DEM 数据 计算 而 
引起 的 精度 损失 。 结 论 如 下 : 

(1) 三 点 拟 合法 计算 的 曲率 值 较 最 小 二 乘法 的 
结果 大 ,空间 分 布 上 存在 明显 差异 ;同一 种 方法 (三 
点 拟 合法 ) 提 取 等 高 线 曲 率 , 选 点 方式 不 同 ,曲率 值 
在 空间 分 布 上 差异 不 明显 。 利 用 三 点 法 计算 的 结 
曲率 值 频 数 分 布 总 体 趋势 一 致 ,并 在 一 定 值 域内 集 
中 ,但 最 小 二 乘法 计算 的 等 高 线 曲 率 值 高 度 集中 于 
(0,0.02] 的 范围 内 。 
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Extraction and analysis of the contour curvature 


based on different algorithms 


WANG Ning', YAO Zhi-hong” 
(1 College of Urban and Environmental Science , Northwest University ,Xi’an 710127 , Shaanxi, China ; 
2 College of Surveying and Geo-Informatics , North China University of Water Resources and Electric Power , 
Zhengzhou 450046 , Henan , China ) 


Abstract; Contour curvature is one of the fundamental topographic indicators that reflects the convexities in the 
horizontal direction of the topographical surface and expresses the divergence and confluence patterns of the motion 
of the surface matters. Calculating the contour curvature can measure the information conversion rate of the DEM in- 
terpolated by the vector contours, and provide a reference for setting the appropriate DEM resolution. This study 
based on vector contour data and DEM of Xiannangou watershed in Ansai County , Shaanxi Province, China, using 
the circle fitting algorithms (adjacent three-point method , interval three-point method and least square method) and 
surface fitting models (Evans model, Shary model and Zevenbergen model) to calculate the contour curvature re- 
spectively. The comparative analysis between the real topography of the site and the computational results indicated 
as follows: (1) The results from the three-point fitting method are different from the results using the least squares 
method, and the former describes the spatial pattern distribution of the contour curvature more accurately. (2) The 
frequency distribution of the results calculated by the least squares method is the highest , and the difference between 
the frequency curves of the two three-point fitting methods is small. (3) In the calculation results of the DEM , the 
results from the S-model display a larger difference than that by the E-model and the Z-model in the spatial pattern. 
The calculation results based on the E-model have a better effect on the contour curvature description. The differ- 
ence presented by this paper when using different algorithms to calculate the contour curvature provides important 
reference in retrieving the contour curvature information in the practical digital terrain analysis. 


Key words: digital terrain analysis; digital elevation model (DEM) ; contour curvature; circle fitting algorithms 


